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基础 研究 


基于 边缘 检测 算 子 的 Huber 正则 化 病 值 选择 方法 在 低 剂量 CT 


重建 中 的 应 用 


张 善 立 , 张 华 , 胡 德 斌 , 曾 栋 , 边 兆 英 , 路 利 军 , 马 建华 ， 


南方 医科 大 学 生物 医学 工程 学 院 , 广 东 广州 510515 


黄 静 


摘要 :目的 研究 两 种 不 同 的 Huber 正 则 化 阔 值 自 适应 选取 方法 及 其 在 低 剂 量 CT 人 迭代 重建 中 的 应 用 。 方 法 针对 低 剂量 CT 重 


建 采用 基于 Huber 正 则 化 的 迭代 重建 技术 ,Huber 正 则 化 阔 值 的 选取 分 别 基 于 全 局 和 局 部 边缘 保持 算 子 。 结 果 仿真 数据 的 实 
验 结果 表明 ,两 类 Huber 正 则 化 阔 值 自 适应 选取 方法 均 能 较 好 地 抑制 重建 图 像 中 的 噪声 和 伪 影 。 结 论 两 类 Huber 正 则 化 阔 值 


自 适 应 选择 方法 均 可 实现 低 剂 量 CT 优质 重建 。 
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Edge-detecting operator-based selection of Huber regularization threshold for low-dose 


computed tomography imaging 


ZHANG Shanli, ZHANG Hua, HU Debin, ZENG Doneg, BIAN Zhaoying, LU Lijun, MA Jianhua, HUANG Jing 
School of Biomedical Engineering, Southern Medical University, Guangzhou, 510515, China 


Abstract: Objective To compare two methods for threshold selection in Huber regularization for low-dose computed 
tomography imaging. Methods Huber regularization-based iterative reconstruction (IR) approach was adopted for low-dose 
CT image reconstruction and the threshold of Huber regularization was selected based on global versus local edge-detecting 
operators. Results The experimental results on the simulation data demonstrated that both of the two threshold selection 
methods in Huber regularization could yield remarkable gains in terms of noise suppression and artifact removal. Conclusion 
Both of the two methods for threshold selection in Huber regularization can yield high-quality images in low-dose CT image 


iterative reconstruction. 


Key words: low-dose CT; iterative reconstruction; Huber regularization; threshold choice 


临床 又 射线 CT 检查 中 的 高 剂量 X 射 线 照射 会 引发 
白血病 ,癌症 和 其 他 多 种 遗传 性 疾病 ”"。 如 何在 不 影响 
临床 诊断 信息 的 条 件 下 ,最 大 限度 地 降低 CT 放射 剂量 
已 成 为 近年 来 CT 领域 吸 竺 解决 的 关键 性 课题 。 当 前 ， 
诸多 策略 相继 被 提出 以 实现 低 剂 量 CT 成 像 癌 ,其 中 降 
低 管 电流 和 (或 者 ) 减 少 曝光 时 间 (mAs) 是 最 为 简单 的 
操作 方法 ”。 然 而 ,该 方法 虽 可 以 一 定 程度 上 降低 辐射 
剂量 ,但 是 将 导致 相应 投影 数据 中 的 量子 噪声 急剧 增 
加 ,表现 为 重建 图 像 中 的 大 量 噪声 和 伪 影 ”。 针 对 此 问 
题 , 近 代 重 建 方 法 (iterative reconstruction, IR ) 被 引入 
至 CT 图 像 重建 中 ,以 实现 优质 的 低 剂 量 成 像 %。 

为 了 有 效 地 消除 噪声 并 保持 图 像 的 结构 和 边缘 信 
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息 , 了 下 网 像 重 建 模型 中 需要 引入 先 验 信息 加 以 约束 , 数 
学 上 此 约束 亦 称 为 正则 化 (也 称 为 惩罚 函数 )。 通 常 , 正 
则 化 多 设计 成 一 个 二 次 惩罚 函数 ,来 惩罚 图 像 某 一 局 部 
区 域 中 不 同 像素 点 间 的 灰 度 差 值 ,其 效果 是 噪声 得 到 抑 
制 而 部 分 边缘 细节 也 会 丢失 名 。 为 了 克服 二 次 函数 的 
上 述 缺 点 ,多 种 边缘 保持 正则 化 被 提出 ” ,其 中 Huber 正 
则 化 最 具有 代表 性 由。Huber 正 则 化 采用 一 个 非 二 次 的 
惩罚 函数 代替 原来 的 二 次 惩罚 函数 , 且 满 足 凸 性 和 对 称 
性 ,使 得 它 能 在 减少 噪声 的 同时 有 效 地 保持 网 像 边 缘 信 
息 。 然 而 ,Huber 正 则 化 的 一 个 主要 缺点 是 国 值 需要 依 
据 具 体重 建 目 标 人 工 选 定 ,大 大 限制 了 其 更 为 广泛 的 实 
际 应 用 。 针 对 此 问题 ,许多 文献 讨论 了 Huber 正 则 化 阔 
值 选取 方法 上 3。 例如 , Yang 等 加 在 Huber 函数 中 引入 
了 另 一 个 国 值 以 便 区 分 图 像 的 边缘 与 伪 影 ;了 Rousseeuw 
等 利用 鲁 棒 性 统计 方法 来 检测 奇异 值 ,再 利用 全 局 边 
缘 检 测算 子 实现 闷 值 选择 ;Black 和 Sapio[“ 利 用 图 像 
的 局 部 统计 特性 ,提出 了 用 局 部 边缘 检测 算 子 进行 阅 值 
选择 ;He 等 定义 了 一 个 比例 因子 ,用 拉 普 拉 斯 卷 积 核 
作为 离散 平滑 算 子 对 效 值 进行 限制 。 鉴 于 上 述 ,本 文 重 
点 研究 基于 全 局 和 局 部 边缘 检测 算 子 的 两 种 不 同 的 
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Huber 正 则 化 阔 值 自 适 应 选取 方法 ,并 将 其 应 用 于 基于 
投影 数据 统计 特征 的 惩罚 加 权 最 小 二 乘 (penalized 
weighted least squares,PWLS) 低 剂量 CT 图像 迭代 重 
建 中 ,其 算法 性 能 通过 数值 仿真 进行 验证 。 


1 方法 
1.1 惩罚 加 权 最 小 二 乘 CT 重建 算法 
在 不 考虑 散射 噪声 的 情况 下 ,CT 的 测量 可 以 表示 
为 一 个 离散 的 线性 系统 : 
y=Hu (1) 
式 中 ,y= (pp,.…, 4s) 为 待 重建 图 像 的 衰减 系 
数 。y=(y, ys.…yw) 表示 测量 得 到 的 弦 图 数据 (经 过 系统 
校正 和 对 数 变换 之 后 的 投影 数据 ), 其 中 ' 代 "代表 矩阵 的 
转 置 , 互 表示 大 小 为 MxN 的 系统 矩阵 ,元 素 瓦 , 代表 像 
素 / 对 投影 射线 的 权 值 。CT 图 像 重建 的 目的 是 从 测量 
数据 y 中 估计 出 衰减 系数 。 
根据 CT 的 测量 模型 (1) 和 最 大 后 验 (maximum a 
posteriori, MAP) 估计 方法 , 带 有 正则 化 项 R(4) 的 
PWLS 图 像 重 建 能 量 方程 可 表示 为 : 
= argmin((—H) 30 Hr) + gR(2) (2) 
其 中 , >， 表示 元 素 为 0 对 角 和 矩阵 ,人 为 弦 图 数据 y 
的 方差 ,98 为 一 个 用 来 控制 正则 化 强度 的 超 参 数 。 
我 们 可 以 利用 多 种 方法 “来 确定 公式 (2) 中 的 参 
数字 ,在 本 文中 ,我 们 使 用 Ma 等 "提出 的 均值 方差 理 


论 : 


07= Texp( 2 员 Top( (= 125) (3) 

式 中 7 为 入射 X 线 强度 , 态 为 在 像素 :处 弦 图 数据 的 
均值 ,52 为 背景 电子 噪声 的 方差 。 
1.2 Huber 正 则 化 

传统 上 , 式 (2) 中 的 正则 化 项 R(z)j 通 过 对 图 像 中 相 
邻 像素 间 灰 度 的 差 值 进行 加 权 求 和 计算 求 得 ,可 以 描 
述 为 如 下 形式 : 

R= RS olk,j)p (wr) (4) 


式 中 ) 为 图 像 域 中 所 有 像素 的 索引 ,S 表 示 二 维 空 
间 第 /个 像素 点 的 邻 域 。 权 重 w(&, 刀 表示 像素 4 和 像素 7 
相互 关系 的 正常 数 ,为 正 值 并 且 是 对 称 的 , 即 w(&, 力 0 
有 目 w(k, 门 =w(j, 有)。w(k, 让 通常 被 设 定 为 邻 域内 像素 
kj 之 间距 离 的 反比 。g 表 示 一 个 凸 的 正 势 函数 并 量 满 
足 g(0) = 二 0。gqg(7) 的 选择 方法 有 许多 种 ,本 文 使 用 Huber 
正则 化 pmwe， 其 定义 为 : 


全 二 22/2; It|<6 
Se alel =052, ds 


(5) 


如 图 1 所 示 ,Huber 正 则 化 为 一 个 凸 函数 , 旦 满足 对 
称 性 , 它 可 以 根据 图 像 当 前 像素 点 与 其 邻 域内 其 他 像素 
点 间 的 灰 度 差 值 来 调整 对 不 连续 区 域 的 惩罚 ,这 一 作用 
由 其 羡 值 6 实现 。 阅 值 5 把 Huber 正则 化 分 为 3 个 部 分 ， 
当 图 像 梯 度 满足 |t|6 时 ,正则 化 的 三 次 项 部 分 对 图 像 
连续 区 域 进 行 较 重 的 惩罚 ; 当 梯 出 直 > 6 时 ,正则 化 的 线 
性 部 分 对 图 像 进行 较 轻 的 惩罚 来 保持 图 像 的 锐利 特性 ， 
比如 图 像 的 边缘 。 

综 上 所 述 ,在 CT 图 像 重建 中 带 有 Huber 正 则 化 的 
PWLS 算 法 (PWLS-Huber) 的 目标 降 数 可 以 写成 以 形式 : 
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图 1 Huber 正 则 化 
Fig.l Example of Huber regularization with 
different thresholds. 


/=arg min(0—H) 30 = BRiwe(p)| (6) 
以 三 
其 中 Huber 正 则 化 项 为 : 
Riwea (1) = Zi wk, )p rwal Ls 一 0 (7) 
j kES; 


在 实际 运用 中 ,我 们 用 一 种 改进 的 共 斩 梯 度 法 来 优 
化 此 目标 函数 。 
1.3 Huber 正 则 化 阐 值 选择 方法 

Huber 正则 化 的 主要 优点 是 它 能 根据 图 像 当 前 像 
素 点 与 其 相 邻 像素 点 间 的 灰 度 差 值 来 调节 对 不 连续 区 
域 的 惩罚 ,这 一 作用 由 其 阔 值 S 实 现 。 因 此 ,对 Huber 郴 
数 来 说 选择 一 个 合理 的 阔 值 是 非常 重要 的 。 针 对 此 问 
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题 ,本 文 主要 研究 基于 全 局 和 局 部 边缘 保持 算 子 的 
Huber 正 则 化 阔 值 选取 方法 以 及 其 在 低 剂 量 CT 图 像 重 
建 中 的 应 用 。 
1.3.1 全 局 边缘 检测 方法 各 向 异性 扩散 被 视 为 一 个 
稳健 统计 学 过 程 , 它 能 从 含有 噪声 的 数据 中 估算 出 
分 段 平 请 的 图 像 ,已 经 在 图 像 边 缘 保持 中 得 到 了 广 
泛 应 用 。 根 据 Perona 和 Malik 的 传统 公式 ,具有 较 
大 梯度 幅 值 的 像素 点 被 视 为 奇异 点 ,又 被 视 为 图 像 
的 边缘 。 这 一 解释 可 用 来 定义 和 检测 图 像 分 段 平 滑 
区 域 的 边界 。 当 一 个 像素 点 的 梯度 幅 值 小 于 阔 值 9 
时 ,我 们 认为 该 像素 点 不 属于 奇异 点 ;否则 ,该 像素 
点 属于 奇异 点 。 全 局 边缘 检测 方法 的 主要 思想 是 阔 
值 6 应 该 描述 整 幅 图 像 数据 的 主要 方差 特性 ”, 即 效 
值 应 该 描述 图 像 中 除 边缘 外 其 他 像素 点 梯度 的 变 
化 。 根 据 这 一 思想 ,我们 用 稳健 统计 学 方式 计算 出 
图 像 所 有 像素 点 梯度 的 方差 ,然后 利用 各 向 异性 扩 
散 方法 确定 空间 变化 的 边缘 阔 值 。 因 此 ,对 于 阔 值 C 
的 鲁 棒 性 测量 可 通过 如 下 方法 : 
66=1.4826MAD(VT) 
= 1.4826median(| VI— mediani(||vI DD 


式 中 ,MAD 表 示 中 位 数 绝对 偏差 。 鉴 于 均值 为 0 

方差 为 1 的 标准 正 态 分 布 的 中 位 数 绝对 偏差 是 0.6745= 
1/1.4826, 故 将 常数 项 设 为 1.4826。 
1.3.2 局 部 边缘 检测 方法 基于 全 局 边缘 检测 算 子 的 
Huber 正 则 化 阔 值 对 于 整 幅 图 像 来 说 是 相同 的 ,然而 我 
们 也 可 以 根据 图 像 的 局 部 统计 特性 来 构造 一 个 可 变化 
的 阔 值 :对 平滑 的 区 域 ,选择 较 大 的 闽 值 以 提高 去 除 噪 
声 和 伪 影 的 能 力 ;对 纹理 较 重 的 区 域 ,选择 较 小 的 阔 值 ， 
以 保持 图 像 的 边缘 。 局 部 边缘 检测 方法 描述 的 是 相对 
于 局 部 图 像 而 言 怎样 通过 图 像 的 梯度 信息 确定 某 一 模 
块 内 的 阔 值 5j( 2 ,gj ,其 主要 思想 是 0 (2 ,9) 应 该 描述 
图 像 中 某 局 部 区 域 数据 的 方差 特性 。61(p,g) 在 整个 图 
像 数 据 中 是 变化 的 , 某 一 点 是 否 属于 图 像 边缘 由 局 部 图 
像 梯度 的 变化 特性 决定 。 

局 部 边缘 检测 方法 和 全 局 边缘 检测 方法 有 相同 的 
稳健 统计 学 理论 基础 。 特 别 地 ,我们 考虑 在 图 像 的 一 个 
大 小 为 n xn 的 模块 内 (这 里 的 模块 是 一 个 图 像 像 素 的 空 
间 位 置 函 数 ,这 种 模块 覆盖 整个 图 像 ) ,计算 其 对 应 的 局 
部 阅 值 61( Pp,9), 公 式 可 以 写 为 : 

01(p,9)=1.4826MAD ,a(VIprisrs) (9) 


(8) 


2 实验 与 分 析 

本 文 使 用 XCAT 体 模 * 生 成 仿真 CT 投影 数据 ,并 
用 上 述 两 种 方法 对 该 投影 数据 分 别 进行 重建 。 实 验 中 ， 
标准 剂量 数据 剂量 为 580 mAs, 低 剂量 为 23.2 mAs。 投 


影 数据 采用 文献 [16] 中 所 述 方法 生成 ,其 中 ,仿真 CT 扫 
描 的 成 像 参数 如 下 : 管 电压 为 120 kVp ,旋转 1 周 采 样 值 
为 1160; 探 测 需 单 元 的 个 数 为 672, 大 小 为 0.8448 mm; 
XX 射线 源 到 探测 器 和 旋转 中 心 的 距离 分 别 为 1361 mm 
和 746 mm; 体 模 大 小 为 512x512; 层 厚 为 3 mm。 实 验 平 
台 为 Intel(R)Pentium(R)CPU G620 @ 2.60 GHz 四 核 
处 理 器 ,3GB 内 存 的 PC 机 ,程序 采用 Matlab 实 现 。 基 于 
本 文 研究 的 两 类 Huber 正 则 化 阔 值 自 适应 选取 方法 的 
PWLS-Huber 算 法 迭代 设 定 为 50 次 , 耗 时 大 约 为 40 min。 

本 文采 用 以 下 两 种 图 像 质量 测度 来 评价 重建 图 像 
质量 :局 部 信 品 比 (local signal-to-noise ratio, ISNR) 和 
相对 均 方 根 误 差 (relative root mean square error, 
rRMSE), 其 计算 公式 分 别 如 下 : 


1 Q 
G2 
ISNR= 2 (10) 
说 起 
pa Hm Hnm 
Q m=1 Q m=1 
Q 2 
Hm 一 KG 
AS (11) 
> HG,m 
m=1 
/表示 低 剂量 投影 数据 重建 出 的 图 像 ,wc 表 示 标 准 


图 像 ,zx 是 感 兴趣 区 域 (Region of Interest ROI) 内 像素 
的 索引 ,Q 是 ROI 内 像素 点 的 数量 。 在 实验 中 ,标准 图 
像 由 在 120 kVp、580 mAs 条 件 下 获得 的 投影 数据 采用 
基于 hann 窗 的 FBP 算 法 重建 得 到 。 

图 2 显示 了 采用 仿真 数据 重建 得 到 的 XCAT 体 模 图 
像 和 低 剂量 CT 图 像 。 图 2A 为 标准 剂量 图 像 ,结构 清晰 ， 
几乎 不 含 任 何 噪声 和 伪 影 ,作为 参照 图 像 。 图 2B 为 在 
23.2 mAs 条 件 下 ( 除 标 准 图 像 外 ,其 余 图 像 的 投影 数据 均 
在 23.2 mAs 条 件 下 获得 ) 获 得 的 投影 数据 采用 FBP 重 建 
得 到 的 图 像 。 和 图 2A 相 比 ,图 2B 中 含有 严重 的 噪声 和 伪 
影 。 图 2C 为 采用 基于 全 局 边缘 检测 算 子 的 PWLS-Huber 
算法 重建 得 到 的 图 像 。 从 图 2C 可 以 看 出 图 像 中 条 形 伪 影 
明显 减少 ,噪声 水 平 大 为 降低 ,图 像 质量 得 到 显著 提高 。 
图 2C 为 采用 基于 局 部 边缘 检测 算 子 的 PWLS-Huber 算 
法 重建 得 到 的 图 像 。 在 图 2D 中 我 们 得 到 了 与 图 2C 相 
似 的 结果 。 由 此 可 以 看 出 ,这 两 种 自 适应 Huber 正 则 化 
贱 值 选 择 方法 , 均 可 得 到 比较 理想 的 结果 一 一 大 部 分 吕 
声 和 伪 影 得 到 了 较 好 的 抑制 。 进 一 步 观察 可 以 发 现 , 虽 
然 图 2C 和 2D 中 伪 影 均 得 到 了 较 好 的 抑制 ,但 是 仍 有 部 
分 伪 影 存在 ,图 像 边缘 的 锐利 度 在 一 定 程度 上 降低 。 相 
比 之 下 ,图 2C 对 图 像 边缘 保持 的 相对 较 好 ,看 起 来 较为 
清晰 ;图 2D 则 对 的 噪声 的 抑制 效果 较 好 。 

为 了 定量 分 析 这 两 种 Huber 正 则 化 阐 值 选择 方法 
的 差异 ,我 们 计算 了 对 应 于 图 2A 中 感 兴趣 区 域 的 ISNR 
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图 2 仿真 投影 数据 重建 出 的 XCAT 体 模 图 像 和 低 剂 量 CT 图 像 

Fig.2 XCAT phantom image and low-dose images reconstructed 
from the simulated projection data. A: Image reconstructed by 
the FBP method with optimal hanning filter from the projection 


© 


data acquired with 580 mAs; B: Image reconstructed by the FBP 
method with ramp filter from the projection data acquired with 
23.2 mAs; C: Image reconstructed by the PWLS-Huber algorithm 
with the global-edge-detecting operator from the projection data 
acquired with 23.2 mAs (8=1.0x10); D: Image reconstructed by 


the PWLS-Huber algorithm with the local-edge-detecting 
operator from the projection data acquired with 23.2 mAs 
(=9x9, 8=1.0x 10). The display option ls [0, 0.03]. 


表 1 对 应 于 图 2A 中 所 示 3 个 ROI 的 ISNR 和 rRMSE 


和 rRMSE, 结 果 如 表 1 所 示 。 可 以 看 出 ,相对 于 FBP 重 
建 图 像 , 采 用 基于 全 局 和 局 部 边缘 检测 算 子 的 
PWLS-Huber 算 法 重建 得 到 的 图 像 ISNR 显著 提高 、 
rRMSE 明 显 减 小 ,具有 更 高 的 增益 。 此 外 ,局 部 边缘 检 
测算 子 在 对 连续 区 域 的 噪声 抑制 方面 的 表现 优 于 全 局 
边缘 检测 算 子 。 

图 3 中 绘 出 了 对 应 于 网 2 中 各 重建 网 像 的 垂直 剖面 
线 ,鉴于 剖面 线 中 含有 512 个 像素 点 , 若 全 部 显示 则 难 
以 区 分 各 重建 方法 , 故 仅 截取 其 中 一 段 ( 即 图 2A 中 红色 
线段 所 示 ) 显 示 , 区 间 为 [330, 360]。 由 图 3 可 以 看 出 ,在 
目标 区 域 基于 两 种 不 同 Huber 正 则 化 阔 值 选择 方法 的 
PWLS-Huber 方 法 的 重建 值 和 标准 值 都 较为 接近 ,但 局 
部 边缘 检测 算 子 比 全 局 边缘 检测 算 子 损失 了 更 多 的 图 
像 分 辨 率 。 

对 于 全 局 边缘 检测 方法 ,Huber 正 则 化 的 浆 值 计算 
出 来 就 固定 不 变 了 。 但 是 ,对 于 局 部 边缘 检测 方法 ， 
个 小 模块 内 的 赋值 是 不 同 的 。 为 了 证 明 模块 大 小 对 重 
建 图 像 质量 的 影响 ,我们 对 局 部 边缘 检测 方法 的 做 了 进 


Tab.1 LSNR and rRMSE measured on the ROIs as indicated by the squares in Fig.2A 


ROI 1 ROI2 ROI3 
Method 
ISNR rRMSE ISNR rRMSE ISNR rRMSE 
FBP 14.78 0.0674 12.90 0.0803 3.59 0.2761 
PWLS-Huber (Global) 64.02 0.0187 66.78 0.0217 17.59 0.0683 
PWLS-Huber (Local) 64.64 0.0186 67.48 0.0216 17.64 0.0680 
ee 比较 明显 ; 随 着 模块 的 增 大 ,不 同 模块 间 闭 值 的 差异 也 
ee TP 逐渐 碱 小 。 图 5 展示 了 在 不 同 模块 大 小 情况 下 ,采用 基 
0.040 1 一 一 PWLS-huber (local) 


es PWLS-huber (global) 
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图 3 对 应 于 图 2A 中 红色 线段 (x=408,y=[330,360]) 的 垂 
直 齐 面 曲线 
Fig.3 Vertical profiles located at the Pixel position x=408 
and y from 330 to 360 as indicated by the red line in Fig.2A. 


一 步 的 实验 。 图 4 展示 了 在 不 同 模块 大 小 情况 下 局 部 
边缘 检测 方法 对 61(p,g) 的 估计 值 。 亮 的 区 域 对 应 图 像 
中 纹理 重 的 区 域 , 具 有 更 高 的 61(p,g) 值 。 可 以 看 出 , 当 
所 选 模块 较 小 时 ,不 连续 区 域 中 不 同 模块 间 阔 值 的 差异 


于 局 部 边缘 检测 算 子 的 PWLS-Huber 算 法 重建 得 到 的 
图 像 。 可 以 看 出 61( 2 ,9) 大 小 变化 对 重建 图 像 的 质量 影 
响 不 明显 。 为 了 进一步 量化 分 析 , 表 2 给 出 了 对 应 网 5 
所 示 图 像 的 不 同 ROI 的 1SNR 和 rRMSE。 可 以 看 出 , 模 
块 大 小 从 5x5 到 17x17 变 化 的 过 程 中 重建 图 像 的 质量 变 
化 很 小 。 当 模块 大 小 为 9x9 时 ,我 们 可 以 得 到 一 幅 质量 


图 4 局 部 边缘 检测 方法 对 应 的 9,(p,g) 值 的 估计 , 亮 的 区 域 对 应 
较 大 的 01( 7， 4) 值 


Fig.4 Local estimate of scale $ z( p, gq) Bright areas correspond to 
larger values of61(p, g). A: n=5x5; B: n=%9; C: n=13x13; D: n=17X17. 
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图 5 采用 基于 局 部 边缘 检测 算 子 的 PWLS-Huber 算 法 重建 得 到 
的 图 像 
Fig.5 Images reconstructed by the PWLS-Huber algorithm with 


the local-edge-detecting method of different patch size. A: n=5x5; 
B: 1=9x9; C: n=13x13; D: n=17x17. The display option is [0, 0.03]. 


表 2 对 应 图 5 所 示 图 像 不 同 ROI 的 ISNR 和 rRMSE 


相对 较 高 的 图 像 。 
3 讨论 与 结论 

针对 Huber 正 则 化 的 阔 值 需要 依据 具体 重建 目标 
人 工 选 定 的 问题 ,本 文 主要 研究 了 基于 全 局 和 局 部 边缘 
保持 算 子 的 Huber 正 则 化 阔 值 选 取 方 法 以 及 其 在 低 剂 
量 CT 图 像 重建 中 的 应 用 ,并 且 对 他 们 做 了 定性 和 定量 
的 测量 与 分 析 。 实 验 结果 表明 局 部 边缘 检测 方法 在 对 
连续 区 域 的 噪声 抑制 方面 的 表现 优 于 全 局 边缘 检测 方 
法 。 这 可 能 是 由 于 由 噪声 引起 的 条 形 伪 影 区 域 和 高 误 
减 区 域 的 86.(p,g) 值 相对 较 大 。 但 另 一 方面 ,从 图 像 过 
度 平滑 的 边缘 可 以 看 出 ,局 部 边缘 检测 方法 比 全 局 边缘 
检测 方法 牺 和 性 了 更 多 的 图 像 分 辩 率 。 综 上 所 述 ,这 两 种 


Tab.2 LSNR and rRMSE measures on the ROIs in Fig.5 with different local patch size 


ROIl ROI2 ROI3 
Patch size 
ISNR rRMSE lISNR rRMSE ISNR rRMSE 
Sx5 64.90 0.0185 67.44 0.0216 17.61 0.0679 
9x9 64.64 0.0186 67.48 0.0216 17.64 0.0680 
13x13 64.32 0.0186 67.39 0.0217 17.63 0.0680 
17x17 64.20 0.0186 67.31 0.0217 17.61 0.0681 


基于 稳健 统计 学 的 Huber 正 则 化 阔 值 选取 方法 对 噪声 
和 伪 影 的 抑制 效果 较 好 ,各 有 优点 ,在 使 用 时 ,可 以 根据 
对 重建 图 像 质量 的 具体 要 求 加 以 选择 。 
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